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Abstract 

We study the critical behavior of period doubling in two coupled one- 
dimensional maps with a single maximum of order z. In particurlar, the effect 
of the maximum-order z on the critical behavior associated with coupling is 
investigated by a renormalization method. There exist three fixed maps of 
the period-doubling renormalization operator. For a fixed map associated 
with the critical behavior at the zero-coupling critical point, relevant eigen- 
values associated with coupling perturbations vary depending on the order z, 
whereas they are independent of z for the other two fixed maps. The renor- 
malization results for the zero-coupling case are also confirmed by a direct 
numerical method. 
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Universal scaling behavior of period doubling has been found in one-dimensional (ID) 
maps with a single maximum of order z {z > 1), 

Xi+i = /(Xi) = 1 - A Z>1. (1) 

For all 2; > 1, the ID map (|T]) exhibits successive period-doubling bifurcations as the non- 
linearity parameter A is increased. The period-doubling bifurcation points A = An{z) 
{n = 0, 1,2, . . .), at which the nth period doubling bifurcation occurs, converge to the ac- 
cumulation point A*{z) on the A axis. The scaling behavior near the critical point A* 
depends on the maximum-order z, i.e., the parameter and orbital scaling factors, 6 and a, 
vary depending on z Therefore the order z determines universality classes. 

Here we study the critical behavior of period doubling in a map T consisting of two 
identical ID maps coupled symmetrically: 

{Xi+i = F{xi, Vi) = f{xi) + g{xi, yi), 
(2) 
Vi+i = F{yi, Xi) = fijji) + giyu x^), 

where f{x) is a ID map (|I]) with a single maximum of even order z {z = 2,4,6,...) at 
X = 0, and g{x, y) is a coupling function. The uncoupled ID map / satisfies a normalization 
condition, /(O) = 1, and the coupling function g obeys a condition, g{x,x) = for any x. 

The quadratic-maximum case {z = 2) was previously studied in Refs. [§-0]. In this 
paper, using the renormalization method developed in Ref . , we extend the results for the 
z = 2 case to all even-order cases and investigate the dependence of the critical behavior 
associated with coupling on the order z. 

The period- doubling renormalization transformation Af for a coupled map T consists of 
squaring (T^) and rescaling {B) operators: 

J\f{T) = BT'^B-\ (3) 

Here the rescaling operator B is: 

(a 0] 

B= , (4) 



a 
2 



because we consider only in-phase orbits (xj = yi for all i). 

Applying the renormalization operator M to the coupled map n times, we obtain the 
n times renormalized map T„ of the form, 



T 



(5) 



Xi+i = Fn{xi,yi) = fn{x) +gn{x,y), 
yi+i = Fnijji.Xi) = fniy) + gniy,x). 

Here /„ and Qn are the uncoupled and coupling parts of the n times renormalized function 
Fn, respectively. They satisfy the following recurrence equations H]: 



X 

fn+l{x) = a/n(/n(-)), 

a 

gn+l{x, y) = a/n(/n(-) + 9n{-, -)) 

a a a 

+a9n{fn{-) + 9n{-, -), /n(-) 

a a a a 

+9n{-, -)) - a/n(/n(-)), 

a a a 



(6) 



(7) 



where the rescaling factor a is chosen to preserve the normalization condition /„+i(0) = 1, 
i.e., a = 1//„(1). Eqs. (||) and (J^) define a renormalization operator TZ of transforming a 
pair of functions, {f,9); 



(8) 



A critical map Tc with the nonlinearity and coupling parameters set to their critical 
values is attracted to a fixed map T* under iterations of the renormalization transformation 
AT, 



(f ) 

Jn+1 


= n 


(f \ 

J n 


^ 9n+l J 




[9n) 



(9) 



Xi+i = F*{xi,yi) = f*{xi) + g*{xi,yi), 
Vt+i = F*{yi,Xi). 

Here {f*,9*) is a fixed point of the renormalization operator TZ with a = which 
satisfies (/*, g*) = 7l{f*, g*). Note that the equation for /* is just the fixed point equation in 
the ID map case. The ID fixed function /* varies depending on the order z Consequently 
only the equation for the coupling fixed function g* is left to be solved. 
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However it is not easy to directly solve the equation for the coupling fixed function. We 
therefore introduce a tractable recurrence equation for a "reduced" coupling function of the 
coupling function g{x,y) 0, defined by 

dg{x,y) 



G{x) = 



dy 



(10) 



y=x 



Differentiating the recurrence equation (0) for g with respect to y and setting y = x, we 
have 



Gn+i{x) = [/;(/„(-)) - 2G„(/„(-))]G„(-) 
a a a 

+G„(/„(-))/;(-). 

a a 



(11) 



Then Eqs. and (|TT]) define a "reduced" renormalization operator TZ of transforming a 
pair of functions (/, G): 



Jn+l 



G. 



n+l 



n 



( f \ 

Jn 

y Gn J 



(12) 



We look for a fixed point {f\G*) of which satisfies {f\G*) = n{f*,G*). Here G* 
is just the reduced coupling fixed function of g* (i.e., G*{x) = dg*{x,y) /dy\y=x)- As in the 
quadratic- maximum case {z = 2) we find three solutions for G*\ 



1. G*(x) 

2. G*{x) 

3. G\x) = ^r\x) - I]. 



0. 



(13) 
(14) 
(15) 



Here the first solution, corresponding to the reduced coupling fixed function of the zero- 
coupling fixed function g*{x,y) = 0, is associated with the critical behavior at the zero- 
coupling critical point, whereas the second and third solutions dependent on the order z are 
associated with the critical behavior at other critical points 0. 

Consider an infinitesimal reduced coupling perturbation (0, <l>(x)) to a fixed point (/*, G*) 
of TZ. We then examine the evolution of a pair of functions, (/*(x), G*{x) + ^(x)) under the 



reduced renormalization transformation TZ. In the linear approximation we obtain a reduced 
linearized operator £, of transforming a reduced coupling perturbation $: 

= [nr(-))-2G*(r(-))]$„(-) 

a a a 

+ [/*'(-) -2G*(-)]<|.„(r(-)). (16) 
a a a 

Here the prime denotes a derivative. If a reduced coupling perturbation $*(a;) satisfies 

z/$*(x) = [C^*]{x), (17) 

then it is called a reduced coupling eigenperturbation with coupling eigenvalue (CE) u. 

We first show that CEs are independent of the order z for the second and third solutions 
(|l^) and ( p!5|) of G*{x). In case of the second solution G*{x) = the reduced linearized 



operator C becomes a null operator, independently of z, because the right hand side of Eq. 
(p!7|) becomes zero. Therefore there exist no relevant CEs. For the third case G*{x) = 
lir'ix) - 1], Eq. becomes 

= $*(-) + $*(/*(-)). (18) 

a a 

When $*(x) is a non-zero constant function, i.e., $*(a;) = c (c : non-zero constant), there 
exists a relevant CE, = 2, independently of z. 

In the zero-coupling case G*{x) = 0, the eigenvalue equation (0) becomes 

u<!>*ix) = /*'(/* (^))r(-) + (/*(-)). (19) 

a a a a 



Relevant CEs of Eq. ([T9|) vary depending on the order z, as will be seen below. 

An eigenfunction ^*(x) can be separated into two components, $*(x) = $*(^^(x) + 

$*(^)(a;) with ^*^^\x) = + ajx + h a*_2a;^~^ and $*(^)(a;) = a*_;^x^~^ + a*a;^ + ■ ■ ■ , 

and the ID fixed function /* is a polynomial in x^, i.e., f*{x) = 1 + c*x^ + Cg^x^^ + ■ ■ ■ . 
Substituting the functions $*, /* and /*' into the eigenvalue equation (|T^, it has the 
structure 



i^al = Y.Mkii{c*})al k,l = 0,1,2,... . (20) 

We note that each (/ = 0, 1, 2, . . .) in the first and second terms in the right hand side of 
Eq. (|19D is involved only in the determination of coefficients of monomials x'^ with k = l + mz 
and k = {z — 1) + mz (m = 0, 1, 2, . . .), respectively. Therefore any a*i with I > z — 1 (in 
the right hand side) cannot be involved in the determination of coefficients of monomials x'^ 
with k < z — 1, which implies that the eigenvalue equation (|20|) is of the form 



^ $*(!) ^ 


^ Ml 




' $*(!) ^ 


^ $*(2) ^ 




M2 ^ 


^ $*(2) ^ 



where Mi is a (z - 1) x (2; - 1) matrix, $*(^) = [a^, a\_^, and $*(^) = (a*_i, a*, . . .). 
From the reducibility of the matrix M into a semi-block form, it follows that to determine 
the eigenvalues of M it is sufficient to solve the eigenvalue problems for the two submatrices 
Ml and M2 independently. 

We first solve the eigenvalue equation of Afi (z/$*(^) = Afi$*(^)), i.e., 

val = J2Mki{{c*})al k,l = 0, . . . , z - 2. (22) 

I 

Note that this submatrix Mi is diagonal. Hence its eigenvalues are just the diagonal ele- 
ments: 

Uk = M,, = ^^ = a'-'-\ k = 0,...,z-2. (23) 

Notice that all u^s are relevant eigenvalues. 

Although z/fc is also an eigenvalue of M, 0) cannot be an eigenvector of M, because 

there exists a third submatrix M3 in M [see Eq. (pi])]. Therefore an eigenfunction $^(a;) 
in Eq. ([T9|) with eigenvalue Uk is a polynomial with a leading monomial of degree k, i.e., 
$*(x) = <^l^^\x) + $fc^^^(x) = alx'' + al_^x'~^ + a*x^ + ■ ■ ■ , where al 0. 

We next solve the eigenvalue equation of M2 (1^$*^^) = M2$*'^^^), i.e., 

= ^MfeK{c*})a;, A;,/ = 2;-l,2;,... . (24) 

I 
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Unlike the case of Mi, (0,$*^^^) can be an eigenvector of M with eigenvalue z/. Then its 
corresponding function $**^^^(a;) is an eigenperturbation with eigenvalue u, which satisfies 
Eq. (p!9|) . One can easily see that ^*^'^\x) = f*'{x) is an eigenfunction with CE = 2, which 
is the 2;th relevant CE in addition to those in Eq. (^). It is also found that there exist an 
infinite number of additional (coordinate change) eigenfunctions = f*\x)[f*^{x) — 

x"] with irrelevant CEs a"" (n = 1,2,...), which are associated with coordinate changes 
||T0| . We conjecture that together with the z (non-coordinate change) relevant CEs, they 
give the whole spectrum of the reduced linearized operator £ of Eq. (p!6|) and the spectrum 
is complete. 

Consider an infinitesimal coupling perturbation g{x,y) [= eip{x,y)] to a critical map at 
the zero-coupling critical point, in which case the map T of Eq. (H) is of the form, 

Xi+i = F{xi,yi) = fA*{x^) +eip{x,y), 
T : i (25) 

1^ y-i+i = F{yi,Xi), 

where the subscript A* of / denotes the critical value of the nonlinearity parameter A and e 
is an infinitesimal coupling parameter. The map T for e = is just the zero-coupling critical 
map consisting of two uncoupled ID critical maps /a*. It is attracted to the zero-coupling 
fixed map (consisting of two ID fixed maps /*) under iterations of the renormalization 
transformation A/" of Eq. (|^). 

The reduced coupling function G{x) of g{x,y) is given by [see Eq. (|10|)] 

d(p{x,y) 



G{x) = e^{x) = e- 



dy 



(26) 



y=x 

Then e^{x) corresponds to an infinitesimal perturbation to the reduced zero-coupling fixed 
function G*{x) = of Eq. (piSf ). The nth image of $ under the reduced linearized operator 
C of Eq. (|T6|) has the form, 

$„(x) = (x) 

2-2 

- + az-i2T'{x) for large n, 

k=0 

(27) 



since the irrelevant part of becomes negligibly small for large n. 

The stability multipliers Ai^„ and A2,n of the 2'^-periodic orbit of the map T of Eq. ( |25D 
are the same as those of the fixed point of the n times renormalized map M^{T) which 
are given by 

^l,n = fni^n), ^2,11 = fn'iS^n) — 2G'„(x„). (28) 

Here (/„, G„) is the nth image of {Ja*, G) under the reduced renormalization transformation 
U [i.e., {fn,Gn) = and x„ is just the fixed point of /„(x) [i.e., x„ = /„(£„)] 

and converges to the fixed point x of the ID fixed map f*{x) as n — oo. In the critical 
case (e = 0), A2,n is equal to Ai^„ and they converge to the ID critical stability multiplier 
A* = f*'{x). Since ^ [C2G]{x) = for infinitesimally small e, A2,n has the form 



V2,n ^ Ai,„ - 2e^n 
rz-2 

^ A* +e 



.k=0 



for large n, (29) 



where = —2ak^l{x) {k = 0, . . . , z — 2) and e^-i = —2az-if*'{x). Therefore the slope 5*, 
of A2,n at the zero-coupling point (e = 0) is 

<9A2,n 



Sn 



z-2 



Y^CkUk + e,_i2'' for large n. (30) 

e=0 k=0 



de 

Here the coefficients {ck ; k = 0, . . . , z — 1} depends on the initial reduced function $(x), 
because a^'s are determined only by $(x). Note that the magnitude of slope Sn increases 
with n unless all e^'s (A; = 0, . . . , z — 1) are zero. 

We choose monomials a;' {I = 0, 1, 2, . . .) as initial reduced functions $(a;), because any 
smooth function $(x) can be represented as a linear combination of monomials by a Taylor 
series. Expressing $( linear combination of eigenfunctions of £2, we have 

$(a;) = = ai^^{x) + a,.if*'{x) 

+Y.Pnr'ix)[rix)-x% (31) 

n=l 

where «/ is non-zero for / < z—1 and zero for / > z—1, and all /?„'s are irrelevant components. 
Therefore the slope Sn for large n becomes 
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eii'V + e^-i2" for I < z - 1, 
Sn^{ ' (32) 
e^_i2" for l>z-l. 

Note that the growth of S'„ for large n is governed by two CEs ui and 2 for I < z — 1 and by 
one CE z/ = 2 for / > ;z - 1. 

We numerically study the quart ic-maximum case {z = 4) in the two coupled ID maps 
(p5|) and confirm the renormalization results (|3^). In this case we follow the periodic orbits of 
period 2" up to level n = 15 and obtain the slopes S'„ of Eq. ( pOD at the zero-coupling critical 
point {A*, 0) {A* = 1.594 901 356 228 . . .) when the reduced function is a monomial 
(/ = 0,1,...). 

The renormalization result implies that the slopes Sn for I > z — 1 obey a one-term 
scaling law asymptotically: 

Sn = d^rl (33) 

We therefore define the growth rate of the slopes as follows: 

ri,n = • (34) 

Then it will converge to a constant ri as n ^ oo. A sequence of ri „ for $(a;) = is 
shown in the second column of Table I. Note that it converges fast to ri = 2. We have also 
studied several other reduced-coupling cases with $(x) = (/ > 3). In all higher-order 
cases studied, the sequences of ri „ also converge fast to ri = 2. 

When / < z — 1, two relevant CEs govern the growth of the slopes Sn- We therefore 
extend the simple one-term scaling law (BBf) to a two-term scaling law: 



Sn = dir^ + d2r2, for large n, (35) 
where |ri| > |r2| . This is a kind of multiple-scaling law [^. Eq. ( |35D gives 

Sn+2 = ^l'S'n+1 — t2Sn, (36) 

where ti = ri + r2 and t2 = rir2- Then ri and r2 are solutions of the following quadratic 
equation, 



r' - hr + t2 = 0. (37) 
To evaluate ri and r2, we first obtain ti and t2 from S^s using Eq. (pG]): 

C C C C C2 C C 

Note that Eqs. p5|)-(|38|) are valid for large n. In fact, the values of tj's and r^'s (z = 1,2) 
depend on the level n. Thus we denote the values of tj's in Eq. ( |38| ) explicitly by tj^„_i's, 
and the values of rj's obtained from Eq. (|37D are also denoted by rj^„_i's. Then each of them 
converges to a constant as n — * oo: 

lim ti^n = U, lim rj,„ = r^, i = 1,2. (39) 



The two-term scaling law (|35D is very well-obeyed. Sequences ri^„ and r2,n for $(x) = 1 
are shown in the third and fourth columns of Table |. They converge fast to ri = 
(a = —1.6903...) and r2,n = 2, respectively. We have also studied two other reduced- 
coupling cases with $(x) = (/ = 1,2). It is found that the sequences ri_„ and r2,„ for 
I = 1(2) converge fast to their limit values ri = a^(2) and r2 = 2(q;), respectively. 
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TABLES 



TABLE I. In the quartic-maximum case (z = 4), a sequence {ri^n} for a one-term scaling law is 

shown in the second column when ^{x) = x^, and two sequences {ri,n} and {r2,n} for a two-term 
scaling law are shown in the third and fourth columns when ^{x) = 1. 





tf\{ fYi\ ,^3 

Vf l^x j — X 


^yX j — i 




n 








r 



i.yyy yzu z 


A QOft AKK Q 

-4.ozy 4iD0 o 


i.yoo 


6 


2.000 009 3 


-4.829422 6 


2.090 


7 


1.999 994 7 


-4.829 409 8 


1.973 


8 


2.000 0004 


-4.829406 8 


2.039 


9 


1.999 999 6 


-4.829405 8 


1.984 


10 


2.000 000 


-4.829 405 5 


2.018 


11 


2.000 000 


-4.829405 5 


1.992 


12 


2.000 000 


-4.8294054 


2.009 
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